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ABSTRACT 

Following on our previous study of an analytic parametric model to describe the baryonic and 
dark matter distributions in clusters of galaxies with spherical symmetry, we perform an SZ 
analysis of a set of simulated clusters and present their mass and pressure profiles. The simu- 
lated clusters span a wide range in mass, 2 x 1O 14 M < M tot (r2oo) < 1 -0 x 10 15 M , and obser- 
vations with the Arcminute Microkelvin Imager (AMI) are simulated through their Sunyaev- 
Zel'dovich (SZ) effect. We assume that the dark matter density follows a Navarro, Frenk and 
White (NFW) profile and that the gas pressure is described by a generalised NFW (GNFW) 
profile. By numerically exploring the probability distributions of the cluster parameters given 
simulated interferometric SZ data in the context of Bayesian methods, we investigate the ca- 
pability of this model and analysis technique to return the simulated clusters input quantities. 
We show that considering the mass and redshift dependency of the cluster halo concentration 
parameter is crucial in obtaining an unbiased cluster mass estimate and hence deriving the 
radial profiles of the enclosed total mass and the gas pressure out to r2oo- 
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1 INTRODUCTION 

Determining the properties of clusters of galaxies such as their to- 
?* tal and baryonic mass offers an independent and powerful cosmo- 
. C-> ' logical tool to constrain the parameters of the ACDM model. The 
' mass distribution of clusters is usually measured using a variety 
5^ , of observational methods, including X-ray, Sunyaev- Zeldovich 
Cp_' (SZ)(Sunyaev & Zeldovich 1970; Birkinshaw 1999; Calrstrom, 
Holder & Reese 2002), and gravitational lensing analyses. These 
methods are often based on parameterised cluster models for the 
distribution of the cluster dark matter and the thermodynamical 
properties of its intra-cluster medium (ICM). However, these var- 
ious approaches usually lead to different estimates of the cluster 
mass. This is due to either extrapolating to halo masses and to red- 
shifts that are not well sampled by the data or fitting for model 
parameters to which the data are insensitive. 

In this letter, we perform a detailed analysis of a sample of 
nine simulated SZ observations of galaxy clusters in a mass range 
of 2 x 10 14 M Q < M tot (r 200 ) < 1.0 x 10 15 M G at redshift z = 0.3 
as if observed with Arcminute Microkelvin Imager (AMI) (AMI 
Consortium: Zwart et al. 2008). To study the cluster total mass, we 
use the model described in Olamaie et al. (2012), which has the 
following main characteristic features: (1) host halo density pro- 
file follows a Navarro, Frenk and White (NFW) (Navarro et al. 
1997) profile and the gas pressure is described by a generalised 
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NFW (GNFW) profile (Nagai et al. 2007) with fixed shape param- 
eters, both in accordance with numerical simulations; (2) the gas 
distribution is in hydrostatic equilibrium with the cluster total grav- 
itational potential dominated by dark matter and both dark matter 
and gas are spherically symmetric; and (3) the local gas fraction is 
much less than unity throughout the cluster, i.e. PgM . , <k 1 for all r. 
This final assumption allows us to write p to t(r) = PdmM +Pws( r ) ~ 
PdmM- We show that assuming the dark matter halo concentra- 
tion parameter as an independent free parameter in the analysis has 
the potential to introduce biases in the cluster mass estimate, and 
hence it is crucial to consider the mass and redshift dependency of 
this parameter in the analysis. Throughout, we assume a ACDM 
cosmology with C1 M = 0.3 , Q.a = 0.7 , cr 8 = 0.8 , h = 0.7 , wo = 
-1, w a = 0. 



2 MODELLING AND ANALYSIS OF SIMULATED 
INTERFEROMETRIC SZ OBSERVATIONS 

As the SZ surface brightness is proportional to the line-of- sight in- 
tegral of the pressure of the hot plasma in the ICM, SZ analysis 
of galaxy clusters provides a direct measurement of the pressure 
distribution of the ICM. 

The observed SZ surface brightness in the direction of electron 
reservoir may be described as 

.dB v 



SI V = T CMB yf(v) 



dT 



(1) 



T=T CMB 
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Here B v is the blackbody spectrum, 7cmb = 2.73 K (Fixsen 
et al. 1996) is the temperature of the CMB radiation, f(v) = 
(jc^tj - 4J (1 + S(x, T c )) is the frequency dependence of thermal 

SZ signal, x = , h p is Planck's constant, v is the frequency 

and k B is Boltzmann's constant. 6(x, T e ) takes into account the 
relativistic corrections due to the relativistic thermal electrons in 
the ICM and is derived by solving the Kompaneets equation up 
to the higher orders (Rephaeli 1995, Itoh et al. 1998, Nozawa 
et al. 1998, Pointecouteau et al. 1998 and Challinor and Lasenby 
1998). It should be noted that at 15 GHz (AMI observing fre- 
quency) x = 0.3 and therefore the relativistic correction, as shown 
by Rephaeli (1995), is negligible for k B T c < 15keV. The dimen- 
sionless parameter y, known as the Comptonization parameter, is 
the integral of the number of collisions multiplied by the mean 
fractional energy change of photons per collision, along the line 
of sight 



y 



o-t p 



n<.(r)k B T e {r)&l ■ 



P e (r)dl, 



(2) 



where n e (r), P c (r) and r c are the electron number density, pressure 
and temperature at radius r respectively, crj is Thomson scattering 
cross-section, m c is the electron mass, c is the speed of light and dl 
is the line element along the line of sight. It should be noted that in 
equation (2) we have used the ideal gas equation of state. 

Moreover, the integral of the Comptonization y parameter over 
the solid angle £1 subtended by the cluster (Fsz) is proportional to 
the volume integral of the gas pressure. It is thus a good estimate 
for the total thermal energy content of the cluster and hence its 
mass (see e.g. Bartlett & Silk 1994). The Ysz parameter in both 
cylindrical and spherical geometries may be described as 



m e c- 



r 



P c {r)2rcs ds, 



P c (r')47ir 2 dr' ', 



(3) 



(4) 



where R is the projected radius of the cluster on the sky. 

In this context we use the model described in Olamaie et al. 
(2012), with its corresponding assumptions on the dynamical state 
of the ICM, to model the SZ signal and determine the radial profiles 
of M 1o1 and P e for nine simulated clusters. The model assumes that 
the the dark matter density follows a Navarro, Frenk and White 
(NFW) profile (Navarro et al. 1997) and the ICM plasma pressure 
is described by the generalised NFW (GNFW) profile (Nagai et al. 
2007), 

Ps 



PdmM 
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where p s is an overall normalisation coefficient, R s is the scale ra- 
dius where the logarithmic slope of the profile d In p(r)/d In r = -2, 
P e i is also an overall normalisation coefficient of the pressure pro- 
file and r p is the scale radius. It is common to define the latter in 
terms of rsoo, the radius at which the mean enclosed density is 500 
times the critical density at the cluster redshift, and the gas con- 
centration parameter, c 500 = r soo /r p . The parameters (a,b,c) de- 
scribe the slopes of the pressure profile at r « r p , r > r p and 
r <K r p respectively. In the simplest case, we follow Araaud et al. 
(2010) and fix the values of the gas concentration parameter and 
the slopes to be (c 50 o,a,b,c) = (1.156, 1.0620,5.4807,0.3292). It 



is also common practice to define the halo concentration parame- 
ter, C200 = -j 2 - . where j'2oo is the radius at which the enclosed mean 
density is 200 times the critical density at the cluster redshift. The 
cluster model parameters: p s , R s and P e i and hence the pressure and 
the integrated mass distributions may be derived under the follow- 
ing assumptions: spherical symmetry; hydrostatic equilibrium; and 
that the local gas fraction is much less than unity, equations (3) to 
(11) in Olamaie et al. (2012). It should be noted that in Olamaie 
et al. (2012), we considered the halo concentration parameter, C200 
as an input free parameter and studied the cluster profiles for a dis- 
tribution of concentrations at a given mass and redshift. However, 
the results of our analysis showed that c 2 oo remains unconstrained. 
C200 is a physical parameter, which reflects the background den- 
sity of the Universe, and so is not a parameter that just defines the 
shape or the slope of the profile that can not be constrained by SZ 
only analysis of galaxy clusters. This therefore suggests that the 
concentration depends on the other physical sampling parameters 
for a given set of cosmological parameters, i.e., mass and redshift. 
Hence, in this letter, we consider such a dependency in our analysis 
and instead derive C2oo- 

We also note that The results of the studies on the relation 
between the structural properties of the dark matter halos such as 
concentration, spin and shape with mass and the redshift from both 
N- body simulations of cosmological structure formation in a Cold 
Dark Matter (CDM) Universe and observations of clusters of galax- 
ies show a clear dependence of the concentration parameter on the 
halo mass and its redshift. This is based on the fact, as first dis- 
cussed by Navarro, Frenk & White (1996) and (1997), that the 
concentration parameter reflects the background density of the Uni- 
verse at its formation time and hence as small objects form first in a 
hierarchical universe, lower mass halos will be more concentrated 
than the massive ones (Navarro, Frenk & White 1997; Eke et al. 
2001; Bullock et al. 2001; Pointecouteau et al. 2005; Maccio et al. 
2007; Vikhlinin et al. 2006; Comerford et al. 2007; Salvador- Sole 
et al. 2007; Buote et al. 2007; Neto et al. 2007; Duffy et al. 2008; 
Mandelbaum et al. 2008; Rudd et al. 2008; Corless et al. 2009; 
Munoz-Cuartas et al. 2011; Ettori et al. 2011; Bhattacharya et al. 
2011). Moreover, all of these studies show that the concentration- 
mass relation is well fitted by a power law over the mass range 
10" - 10 15 /7~'M Q . However, studies on determining the normalisa- 
tion coefficient and the slope of such relation are still ongoing, and 
in particular, cluster observations have yet to investigate this fur- 
ther. In this letter, we, therefore, decided to use the relation derived 
by Neto et al. (2007) from W-body simulation which has also been 
used by Corless et al. (2009) in their weak lensing analysis of three 
real clusters, namely 

, -0.1 



c'200 



5.26 / M lot (r 20 o) 



1 +M 10 14 /i-'M o 



(7) 



We also studied the c\ n - M ylI relation given by Munoz- Cuartas 
et al. (201 1) where the normalization coefficient and the slope vary 
with cosmic time. However, the results were the same as using 
equation (7) within our cluster halo mass range. 

We generate a sample of nine simulated SZ clusters equally 
spaced in the mass range 2 x 10 14 M o < M tot (r 2 oo) < 1.0 x 10 15 M Q 
using the above mentioned model, equation (7) and the input pa- 
rameters M lot (r 2 oo), z and /g as (r 2 oo) listed in Tab. 1; this set of pa- 
rameters fully describes the Comptonization y parameter. Further 
details of generating simulated SZ skies and observing them with a 
model AMI small array (SA) are described in Hobson & Maisinger 
(2002), Grainge et al. (2002), Feroz et al. (2009) and AMI Consor- 
tium: Olamaie et al. (2012). 
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Table 1. M tol (r2oo) used to generate simulated clusters and the thermal noise 
levels reached in the simulated observations of the cluster sample. Here <tsa 
refers to the thermal noise levels reached in SA maps. All the clusters are 
generated at fixed redshift z = 0.3 and fixed / g (r2oo) = 0.13. . 
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Figure 1. ID marginalised posterior distributions of sampling parameters of 
clsim2 when C200 is also assumed to be an input parameter. Green vertical 
lines are the true cluster parameter values as given in Tab. 1 and the red 
vertical lines are the mean of the probability distributions of the parameters. 

The sampling parameters in our Bayesian analysis are C = 
(x c ,y c , M tot (r2oo), /g(/2oo), z), where x c and y c are cluster projected 
position on the sky. We further assume that the priors on sampling 
parameters are separable (Feroz et al. 2009) such that 

7r(0 c ) = tt(x c ) n(y c ) K{M T {r 20 o)) n(f g (r 2 oo)) n(z). (8) 

We use Gaussian priors on cluster position parameters, centred on 
the pointing centre and with standard-deviation of 1 arcmin and 
adopt a S function prior on redshift z- The prior on M tot (r 2 oo) is 
taken to be uniform in logM in the range M rain = 10 14 M Q to M max = 
6x 10 15 M G and the prior of / gas (r20o) is set to be a Gaussian centred 
at the / gas = 0.13 with a width of 0.02 (Vikhlinin et al. 2005, 2006; 
Komatsu et al. 2011; Larson et al. 2011). It should be noted that 
for the two low mass clusters we set the minimum mass to M ra i„ = 
0.4 x 10 14 M Q in the prior range. A summary of the priors and their 
ranges are presented in Tab. 2. 

Further details of our Bayesian methodology, modelling in- 
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Figure 2. ID marginalised posterior distributions of sampling parameters 
of clsim2 when C200 is calculated using equation (7). Green vertical lines 
are the true cluster parameter values as given in Tab. 1 and the red vertical 
lines are the mean of the probability distributions of the parameters. 



Table 2. Summary of the priors on the sampling parameters. Note that 
N(ji, o") represents a Gaussian probability distribution with mean /i and stan- 
dard deviation of cr and U(a, b) represents a uniform distribution between a 
and b. 



Parameter 


Prior 


Xc , >'c 


N(0 , 60)" 


log M to ,(r 2 oo) 


U(U , 15.8) M Q 


/gasC^Oo) 


#(0.13 , 0.02) 



terferometric SZ data, primordial CMB anisotropies, and resolved 
and unresolved radio point-source models are given in Hobson & 
Maisinger (2002), Feroz & Hobson (2008) and Feroz et al. (2009), 
AMI Consortium: Davies et al. (2011) and AMI Consortium: Ola- 
maie et al. (2012). 



3 RESULTS AND DISCUSSION 

Fig. 1 shows ID marginalised posterior distributions of sampling 
parameters for clsim2 (one of the lowest mass clusters) when C200 
is also assumed to be a sampling parameter, (Olamaie et al. 2012). 
The green vertical lines are the true values of the simulated cluster 
parameters while the red vertical lines represent the mean values of 
the probability distributions. It should be pointed out that similar 
results are obtained for all clusters in our sample. This results show 
that while this form of the parameterization can constrain cluster 
projected position on the sky and M 10( (r 2 oo), the halo concentration 
parameter, c 2 oo, is unconstrained and its mean value is just the mean 
of the prior range which suggests that the mean value of the distri- 
bution is strongly driven by the prior. This result, as we have seen 
in our previous studies of clusters of galaxies (AMI Consortium: 
Olamaie et al. 2012), also suggests a correlation between C200 and 
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Figure 3. Integrated mass (left) and pressure (right) profiles as a function of r for nine simulated SZ clusters. In each panel, the plots with * show the profiles 
when C200 is assumed as an input (sampling) parameter and the plots with o show the profiles when c 2 oo is calculated using equation (7). 



the other two physical sampling parameters, i.e. M lot (r 20 o) and z. 
The existence of such a correlation means that the analysis may re- 
sult in a biased estimate of cluster parameters if it is not considered 
in the analysis. 

Fig. 2 shows ID marginalised posterior distributions of sam- 
pling parameters for clsim2 when we take into account the depen- 
dency of halo concentration on both the formation time and the 
dynamical state of the halo using equation (7). From the results it 



is clear that this form of the parameterization can constrain c 2 oo as 
well as other cluster parameters. Moreover, while NFW profiles are 
usually fitted using the two parameters R s and c 2 oo this parametrisa- 
tion makes the profile a one-parameter profile. We notice that simi- 
lar results were obtained upon analysing all the clusters in our sam- 
ple. We also note that / g (r 2 o ) is hardly constrained in both param- 
eterisations indicating that the gas fraction can not be constrained 
using SZ only data. 
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The left panel in fig. 3 presents the integrated mass profiles for 
our sample of nine SZ simulated galaxy clusters and the right panel 
shows the pressure profiles of these clusters. In each panel we have 
plotted the profiles out to r 20 o using the two forms of parameteri- 
zations: (1) assuming C200 as an input parameter in the analysis (*) 
and (2) calculating c 2 oo using equation (7), (o). 

From the plots, the difference in estimating the cluster mass 
and its gas pressure using two forms of parameterizations is clear 
and becomes more significant as M tot (r 20 o) increases. Our result- 
ing constraints on c 2 oo indicate that because this parameter is com- 
pletely unconstrained when assumed as a sampling parameter, its 
best fit value is always the mean of the assumed prior range. This 
strong dependency on the prior range may indeed lead to a biased 
estimate of the cluster parameters including its mass. 



4 CONCLUSION 

We have studied the recovery of M tol (r) and P c {r) from the SZ effect 
for a sample of nine simulated galaxy clusters (2.0 x 10 14 M o < 
M tot (r 2 oo) < 1.0 x 10 15 M G ) using the model described in Olamaie 
et al. (2012). This is motivated by the fact that SZ surface brightness 
is proportional to the line of sight integral of the ICM plasma so that 
SZ data can potentially constrain the cluster total mass. 

To obtain an unbiased mass estimate we have carried out a de- 
tailed analysis of a series of simulated clusters using two different 
parameterizations within our model and its corresponding assump- 
tions (Olamaie et al. 2012). In the first parameterization we assume 
that the halo concentration parameter c 2 oo is also a sampling param- 
eter. However, the results of the analysis show that the simulated SZ 
data can not constrain this parameter and therefore its mean value 
is driven by the prior range which may introduce biases in the ul- 
timate cluster mass estimate. This results also suggests that c 2 oo 
depends on the other model parameters. 

In the second parameterization we consider the correlation of 
c 2 oo with M tot (r 2 oo) and z within ACDM Universe (equation 7) as 
higher mass halos that are forming today are less concentrated than 
halos of lower mass that built up at an earlier epoch, where the 
mean density was higher. This parameterization clearly constrains 
c 20 o as AMI SZ data can constrain M tot (r 20 o). We hence conclude 
that in order to obtain a robust estimate on cluster physical parame- 
ters including its mass it is crucial to consider the mass and redshift 
dependency of c 20 o as precisely as possible. 



ACKNOWLEDGMENTS 

The authors thank their colleagues in the AMI Consortium for nu- 
merous illuminating discussions regarding the modelling of galaxy 
clusters. The data analyses were carried out on the COSMOS UK 
National Supercomputer at DAMTP, University of Cambridge and 
we are grateful to Andrey Kaliazin for his computing assistance. 
MO acknowledges an STFC studentship. 



Bhattacharya S., Habib S„ Heitmann K, 2011, arXiv: 11 12.5479 
Birkinshaw M., 1999, PhR, 310, 97 

Bullock J. S., Kolatt T. S., Sigad Y., Somerville R. S., Kravtsov 
A. V., Klypin A. A., Primack J. R., Dekel A., 2001, MNRAS, 
321,559 

Buote D. A., Gastaldello F., Humphrey P. J., Zappacosta L., Bul- 
lock J. S., Brighenti F., Mathews W. G., 2007, ApJ, 664, 123 
Carlstrom J. E., Holder G. P., Reese E. D., 2002, ARA &A, 40, 643 
Challinor A., Lasenby A., 1998, ApJ, 499, 1 
Comerford J. M., Natarajan P., 2007, MNRAS, 379, 190 
Corless V. L., King L. J., Clowe D„ 2009, MNRAS, 393, 1235 
Duffy A. R., Schaye J., Kay S. T, Dalla Vecchia C, 2008, MNRAS, 
390, L64 

Eke V. R., Navarro J. E, Steinmetz M., 2001, ApJ, 554, 114 
Ettori S., Gastaldello E, Leccardi A., Molendi S., Rossetti M., 

Buote D., Meneghetti M., 2011, A&A, 526, 1 
Feroz E, Hobson M. P., 2008, MNRAS, 384, 449 
Feroz E, Hobson M. P., Bridges M., 2009, MNRAS, 398, 1601 
Feroz E, Hobson M. P., Zwart J. T. L., Saunders R. D. E., Grainge 

K. J. B., 2009, MNRAS,398, 2049 
Fixsen D. J., Cheng E. S., Gales J. M., Mather J. C, Shafer R. A., 

Wright E. L., 1996, ApJ, 473, 576 
Grainge K., Jones M. E., Pooley G, Saunders R., Edge A., Grainger 

W. E, Kneissl R., 2002, MNRAS, 333, 318 
Hobson M. P., Maisinger K, 2002, MNRAS, 334, 569 
Itoh N., Kohyama Y„ Nozawa S„ 1998, ApJ, 502, 7 
Komatsu E., et al., 2011, ApJS, 192, 18 
Larson D., et al., 2011, ApJS, 192, 16 

Maccio A. V., Dutton A. A., van den Bosch F. C, Moore B., Potter 

D., Stadel J., 2007, MNRAS, 378, 55 
Mandelbaum R., Seljak U, Hirata C. M., 2008, JCAP, 8, 6 
Munoz-Cuartas J. C, Maccio A. V., Gottlober S., Dutton A. A., 

2011, MNRAS, 411,584 
Nagai D., Kravtsov A. V., Vikhlinin A., 2007, ApJ, 668, 1 
Navarro J. E, Frenk C. S., White S. D. M., 1996, ApJ, 462, 563 
Navarro J. E, Frenk C. S., White S. D. M., 1997, ApJ, 490, 493 
Neto A. E, et al., 2007, MNRAS, 381, 1450 
Nozawa S., Itoh N., Kohyama Y., 1998, ApJ, 508, 17 
Olamaie M., Hobson M. P., Grainge K. J. B., 2011, 

arXiv: 1109.2834 
Pointecouteau E., Giard M., Barret D., 1998, A&A, 336, 44 
Pointecouteau E., Arnaud M., Pratt G. W., 2005, A&A, 435, 1 
Rephaeli Y, 1995, ARA &A, 33, 541 

Rudd D. H, Zentner A. R., Kravtsov A. V., 2008, ApJ, 672, 19 
Salvador-Sole E., Manrique A., Gonzalez-Casado G, Hansen S. H, 

2007, ApJ, 666, 181 
Sunyaev R. A., Zeldovich Y. B„ 1970, CoASP, 2, 66 
Vikhlinin A., Markevitch M., Murray S. S., Jones C, Forman W., 

Van Speybroeck L., 2005, ApJ, 628, 655 
Vikhlinin A., Kravtsov A., Forman W., Jones C, Markevitch M., 

Murray S. S., Van Speybroeck L., 2006, ApJ, 640,691 



REFERENCES 

AMI Consortium: Davies et al., 2011, MNRAS, 415, 2708 
AMI Consortium: Olamaie M. et al., 2012, MNRAS, 419, 2921 
AMI Consortium: Zwart J. T. L., et al., 2008, MNRAS, 391, 1545 
Arnaud M., Pratt G. W., Piffaretti R., B ohringer H, Croston J. H„ 

Pointecouteau E., 2010, A&A, 517, A92 
Bartlett J. G, Silk J., 1994, ApJ, 423, 12 



